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Abstract. Topological subsystem codes proposed recently by Bombin are quantum 
error correcting codes defined on a two-dimensional grid of qubits that permit reliable 
quantum information storage with a constant error threshold. These codes require 
only the measurement of two-qubit nearest-neighbor operators for error correction. In 
this paper we demonstrate that topological subsystem codes (TSCs) can be viewed 
as generalizations of Kitaev's honeycomb model to 3-valent hypergraphs. This new 
connection provides a systematic way of constructing TSCs and analyzing their 
properties. We also derive a necessary and sufficient condition under which a syndrome 
measurement in a subsystem code can be reduced to measurements of the gauge 
group generators. Furthermore, we propose and implement some candidate decoding 
algorithms for one particular TSC assuming perfect error correction. Our Monte Carlo 
simulations indicate that this code, which we call the five-squares code, has a threshold 
against depolarizing noise of at least 2%. 
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1. Introduction 

Active quantum error correction [1], Ej [3l H] protects quantum information against noise 
by performing measurement to detect errors, and by correcting them. The goal of active 
error correction is to be able to store and manipulate quantum information arbitrarily 
long in the presence of noise provided that the noise level is below some critical threshold. 

The stabilizer formalism [5] provides a convenient framework for studying such error 
correcting codes. A special subclass of these codes are topological stabilizer codes P [7] . 
Topological stabilizer codes, in particular the surface or toric code family [6l |8] offer 
practical advantages over other quantum error correcting codes. First, qubits can be 
laid out on a two-dimensional grid and error correction measurements can be performed 
in-place on the grid. Secondly, topological stabilizer and topological subsystem codes 
[9] have noise thresholds which correspond to phase transitions in associated statistical 
mechanics models; these thresholds asymptote to a finite value in the (thermodynamic) 
limit of large block size, obviating the need for code concatenation. Thirdly, topological 
stabilizer and subsystem codes may allow for the implementation of various (but not all) 
gates such as the important CNOT gate by topological 'code deformation' |TOl [TTl [T2]. 
This topological implementation of logical gates does not require an additional overhead 
in the number or layout of qubits. Universality is obtained by supplementing the 
topological implementation of logical gates by magic state distillation [13] for those gates 
which cannot be implemented topologically. An important advantage of this topological 
implementation of gates is that the tolerable noise rate for the code architecture will 
be largely determined by the quantum memory noise threshold. The quantum memory 
noise threshold for the surface code is among the highest of any codes that has been 
studied [H]. 

In this paper, we will refer to and calculate thresholds for the depolarizing error 
model. In this model X, Y, and Z errors occur on any given qubit with probability 
p/3, and no error occurs with probability 1 — p. The depolarizing probability is p. For 
noisy error correction on the surface code, the quantum memory noise threshold which 
has been achieved using explicit decoding is 0.75 — 1.1% (see e.g. [151 El [161 [T7] . 

In addition, we will only consider noise-free error correction. This means that we 
assume that the syndrome measurements, -obtained for topological subsystem codes 
via the measurement of gauge group generators-, are noise free. With noise-free error 
correction the noise threshold of the surface code against depolarizing noise can be as 
high as 15 — 16% [TH [191 [16]. This number is close to the theoretical maximum in |20j . 

Although the surface code is among the best code candidates developed so far, 
the actual experimental implementation of a surface code architecture is an extremely 
daunting task. Currently, Josephson junction qubits whose coherence times and 
capabilities have progressed rapidly over the past few years, seem to be well-suited 
for a 2D surface code architecture [21]. With an anticipated elementary gate time of 
r = (9(10^®) seconds for Josephson junction qubits, a T2 decoherence time of O(10~^) 
seconds would lead to an effective noise rate of T/T2 = 0(10^^) below the noise threshold 
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of the surface code. 

Beyond a sufficiently-low noise rate, a successful surface code architecture would 
have other technological requirements such the existence of a relatively fast measurement 
and the processing of syndrome data. Since the stabilizer group of the surface code 
is generated by four-body operators, syndrome measurement has to be performed by 
coupling a single ancilla qubit to four neighboring system qubits. In a Josephson junction 
architecture, this could for example mean that a single Josephson junction qubit would 
have to couple to four different resonators. The goal of this work is to ask whether one 
can simplify the physical implementation of error correction by employing topological 
subsystem codes [HI ES] in which 2-body (instead of 4-body) operators need to be 
measured for error correction. The formalism of subsystem codes and their relation 
to stabilizer codes were introduced and analyzed in |23l |2^. Some subsystem codes 
such as the Bacon-Shor code family [211 EH] are characterized by a stabilizer group 
which has non-local generators. These codes, although interesting, are not believed to 
have an asymptotically non- vanishing threshold, unlike the topological subsystem codes 
introduced by Bombin [9]. 

After briefly introducing some stabilizer and subsystem code notation in Section [21 
Section [3] discusses ways of constructing subsystem codes with interesting properties. 
This method includes Bombin's topological subsystem codes, but can also give rise to 
new families of codes. Section 13.11 lists the desired code properties. Section 13.21 sets 
the stage for our results by considering a 'trivial' topological subsystem code based on 
Kitaev's honeycomb model (this code has no logical qubits). This code serves as a 
template for all code constructions introduced later. It also provides us with insights 
in how and under what assumptions the syndrome measurement can be implemented 
by measuring 2-body gauge group generators, see Section [3^ and Appendix A. We will 
define extensions of the honeycomb model to 3-valent hypergraphs and explain how these 
extensions can be used to derive topological subsystem codes, see Sections 13.31 and 13.41 
Our general construction is illustrated with two specific examples in Section 13.51 where 
we define the square-octagon code and the five-squares code. Finally, in Section H] 
we describe a possible decoding algorithm for the five-squares code. We evaluate the 
performance of the decoding algorithm in Section |5l 

2. Background 

Quantum error correcting codes encode quantum information into some subspace C, 
-the code space-, of a larger Hilbert space. After the encoded information is subjected 
to noise, a syndrome measurement is performed to diagnose the error that occurred. 
Finally, the syndrome information is used to return the code into a state that corresponds 
to the originally encoded information. 
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2.1. Stabilizer Codes 

We consider a system comprised of n qubits. Let Xj, Yj, and Zj represent the Pauli 
operators acting on the j-th qubit. Let Vn = {i I, Xi, Zi, X„, Z„) be the Pauh group 
on n qubits. A stabihzer group S is an Abehan subgroup of Vn which does not contain 
— /. One can always choose independent generators S'l, 5*2, . . . Sn-k G S for some k < n. 

The codespace C is the 2^^ dimensional space stabilized by S, i.e., the subspace 
of states for which Sj\ip) = {ip) holds for all j = 1, . . . ,n — k. Such stabilizer 
code encodes k logical qubits whose logical operators are denoted as {Xj, ... fc. 

The centralizer C{S) in Vn includes all Pauli operators in Vn that commute with every 
element of S. For any stabilizer code one has the identity C(S) = {S, Xi, Zi, . . . , Xk, Zk). 
Hence the logical operators of the code are elements of C{S) \S. 

2.2. Subsystem Codes 

Subsystem codes are stabilizer codes in which some logical qubits do not encode any 
information. The presence of these unused logical qubits a.k.a. gauge qubits often allows 
one to simplify encoding and decoding circuits since the initial and final state of the 
gauge qubits is irrelevant. A subsystem code derived from a stabilizer code S is fully 
characterized by its gauge group Q that includes all stabilizers S and all logical operators 
of S acting only on the gauge qubits. The gauge group uniquely defines the stabihzer 
group (up to overall phase factors) such that {iI,S) = Q H C{Q). We shall often ignore 
the overall phase factors and use a simpler identity S = GnC^Q). A subsystem code has 
two species of logical operators. Bare logical operators are the logical operators of S 
that act trivially on the gauge qubits. Multiplying bare logical operators by elements of 
Q one obtains dressed logical operators. Non-trivial bare and dressed logical operators 
are elements of C{Q)\Q and C{S)\Q respectively. Errors are detected and corrected 
by measuring the gauge group generators Gi and deducing the ±1 eigenvalues of the 
generators of S (the syndrome bits) from these measurements. These syndrome bits are 
then used to correct errors modulo operators in the gauge group Q (which do not affect 
the encoded information), see Section [3^ and Appendix A for more details. 

3. Topological Subsystem Codes 

3.1. Desired Code Properties 

The purpose of this section is to introduce a toolbox that allows one to design new 
2D topological subsystem codes. The physical qubits of such codes occupy sites of a 
two-dimensional lattice, whereas the gauge group is generated by spatially local two- 
qubit generators. We begin by introducing some notation and describing desired code 
properties. Let £ = (V, E) be the underlying lattice embedded into a plane or a torus 
with a set of sites V and a set of links (edges) E. Each site u E V represents a physical 
qubit. Pauli operators acting on u will be denoted X„, Yu, and Z^. Each link [u, v) E E 
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represents a gauge group generator Ku^y G Vn acting on a pair of qubits u, v. We shall 
refer to the generators K^^^ as link operators. For example, XuX^, YyYy, and XyYy 
are admissible link operators. For simplicity let us assume that the link operators are 
hermitian, so that = I for any link e & E. Link operators generate a gauge group 

g={K, eeE), 

which in turn specifies a subsystem code as described in Section 12.21 Let us say that a 
Pauli operator P G "Pn is spatially local iff the support of P is confined to a disk of radius 
0(1). Our goal is to construct subsystem codes combining the following properties: 

(CI) The stabilizer group S = Q H C{Q) has spatially local generators. Elements of S 
can be identified with closed homologically trivial loops on the lattice. 

(C2) Syndrome extraction can be implemented by measuring eigenvalues of the two-qubit 
link operators. 

(C3) The code encodes one or more logical qubits. Bare logical operators (in C{Q)\Q) 
can be identified with closed homologically non-trivial loops. 

Properties (CI) and (C3) are often used to describe topological quantum order in the 
framework of stabilizer codes, see for instance [26]. Property (C2) is the main advantage 
of subsystem codes compared to topological stabilizer codes such as the surface code. 
It says that the syndrome extraction requires only two-qubit measurements applied to 
nearest-neighbor sites of the lattice. In contrast, the syndrome extraction for all known 
topological stabilizer codes requires at least four-qubit measurements (it seems hard if 
not impossible to construct a nontrivial stabilizer code with 3- or 2-body, spatially local, 
commuting generators). We shall refer to subsystem codes satisfying conditions (Cl-3) 
as topological subsystem codes. 



3.2. Kitaev's Honeycomb Model on 3-valent Graphs 

Our starting point is Kitaev's honeycomb model [27j and its generalization to arbitrary 
3-valent graphs. It will provide us with a supply of 'trivial' subsystem codes that satisfy 
conditions (CI) and (C2) but do not encode any logical qubits. These trivial codes will 
serve as a template for more sophisticated code constructions introduced later. 
Recall that the honeycomb model [27] is described by a Hamiltonian 

H Jx ^ ^ X,f^X,u Jy ^ ^ YyYy 

{u,v)gEx {u,v)€Ey 

~ Jz ^ Z^Zy (1) 

where qubits occupy sites of the honeycomb lattice C = (V, E) embedded into a torus 
with three species of edges, E = E^UEyUEz, as illustrated in Figure [H We shall regard 
the terms in the Hamiltonian Eq. ([l]) as link operators of a subsystem code, that is, 

{XuXy if (u, v) G Ex, 
Y^Y, if iu,v)eEy, (2) 
ZuZy if {u,v) G E^. 
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Figure 1. Three types of link operators in the Kitaev's honeyconib modeL This 
definition can be extended to the entire lattice using the translation invariance. 



Note that the hnk operators are not independent. Since each site has exactly three 
incident links contributing Pauli operators X, Y, and Z acting on this site, one has an 
identity 

n ~ /. (3) 

e€E 

Furthermore, a pair of link operators K^, K^' anticommutes iff the links e and e' share 
exactly one end-point and commute in all remaining cases, that is, 

KeKe' = {-ly'^^'^'^ Ke'Ke for all e, e' G E, (4) 



where 



r]{e,e) 



1 if e and e' share exactly one site, 
otherwise. 



Consider now a loop operator W{j) = Hee-y associated with any closed loop 7 C 
on the lattice. Using Eq. (j4]) one can easily check that ^^(7) commutes with all link 
operators. Thus, if we consider a subsystem code with a gauge group Q = {K^, e E E), 
then the loop operators are stabihzers of Q, that is, W{'y) E S = Q nC(^) for any closed 
loop 7. In fact, as we prove below (see Lemma [3.11) the stabilizer group S is generated 
by the loop operators. It shows that the honeycomb model provides a natural supply 
of codes satisfying condition (CI). 

Let us now discuss how the syndrome measurement can be implemented by 
measuring eigenvalues of individual link operators. Consider first the extraction of a 
single syndrome bit corresponding to some stabilizer S E S. Suppose we can represent 
S as an ordered product of link operators 

S = Krn---K2Ki where Kj = Ke,, (5) 

such that each Kj commutes with the product of all link operators that appear before 
it, that is, 

[Kj, Kj_i . . . ifi] = for all j = 2, . . . , m. (6) 

Note that Eq. ([6]) implies that S* is a hermitian operator since the total number of 
anticommuting pairs among Ki, . . . ,Km is even. In order to measure the eigenvalue 
a of S, one can perform the nondemolition eigenvalue measurement of Ki, . . . , in 
the order from the left to the right, obtaining outcomes 71, ... ,7^ G {+1, — !}• The 
deduced eigenvalue of 5* is then computed as a = 71 ■ ■ ■ 7m- Let us check that the 
above procedure is equivalent to the nondemolition eigenvalue measurement of S modulo 
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gauge operators. For example, if m = 2 then Eq. implies that Ki and K2 mutually 
commute. Accordingly, the eigenvalues of Ki and K2 can be measured simultaneously. 
Since S = K2K1, the eigenvalue of S is determined by a = 7172. We can now use 
induction to analyze the case m > 2. Define a truncated version of the stabilizer S, 
namely, Sj = Kj ■ ■ ■ Ki. Our induction hypothesis is that the reduced state obtained 
immediately after the measurement of Kj is an eigenvector of 5*^- with an eigenvalue 
(Tj = 7i---7j. Indeed, we have already checked this for j = 1,2. Since Eq. 
implies that Kj^i commutes with Sj, the eigenvalue measurement of -ft'j+i preserves 
the eigenvalue of Sj. Hence the eigenvalue of Sj+i after the measurement of -ft'j+i 
is cTj+i = 7j+icrj. In Appendix A we prove that the existence of the decomposition 
Eq. ([5]) satisfying Eq. (Q is in fact a necessary and sufficient condition under which the 
eigenvalue measurement of S can be simulated by the measurements of the gauge group 
generators. This condition applies to any subsystem code and any fixed choice of the 
gauge group generators. 

Consider now a loop 'j E associated with the boundary of some hexagon. 
Let 61,62,..., 66 be the six consecutive links comprising 7, and Kj = K^^ be the 
corresponding link operators. In order to measure eigenvalue of the loop operator 
14^(7), one can use a decomposition W{'y) = KqK4^K2K^K^Ki, which clearly satisfies 
Eq. dnj. In fact, since the Kj pairwise commute for all odd (even) j, one can measure 
the eigenvalue of W^(7) in only two rounds: Ki, K3, K^ are measured in the first round, 
and K2,K4^,Kq are measured in the second round. Let us point out that the above 
procedure breaks down for loops 7 with an odd length. In fact, one can easily check that 
condition Eq. ([6]) cannot be satisfied for odd loops, and thus the corresponding loop 
operators ^(7) cannot be measured using only two-qubit measurements. It shows that 
bipartiteness of the lattice £ plays a crucial role since for a bipartite graph any closed 
loop has even length. 

Syndrome measurements associated with a pair of adjacent hexagons may interfere 
with each other since one link operator is shared by both stabilizers. However one can 
arrange all the hexagons in 0(1) layers such that hexagons from the same layer do not 
have common links. Now the syndrome extractions for different hexagons in the same 
layer do not interfere with each other and one needs only two measurement rounds 
per each layer. The full syndrome extraction thus requires 0(1) rounds of two-qubit 
measurements 0. It shows that the honeycomb model provides a natural supply of codes 
satisfying conditions (Cl-2). This discussion does not address completely how efficient 
the measurement procedure is from a practical point of view, i.e. precisely how many 
rounds of measurements are needed. 

Unfortunately, as we demonstrate below, the honeycomb model does not satisfy 
condition (C3), that is, the corresponding subsystem code has no logical qubits. To 
set the stage for our results, it will be convenient to analyze a generalization of the 

I For simplicity we ignored the non-local stabilizers associated with homologically non-trivial loops on 
the lattice. Using the same arguments as above, one can measure the eigenvalue of any such stabilizer 
using only two rounds of two-qubit measurements. 
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honeycomb model Eqs. (fT|2]) defined on an arbitrary 3-valent graph C = {V,E). As 
before, we assume that the physical qubits live at sites u & V and the link operators 
Ke obey the commutation rules Eq. (jlj). We note that the commutation rules Eq. (jl]) 
uniquely specify the hnk operators up to local unitaries. Indeed, let u E V be any site 
and f,g,hEEhe the three links incident to u. From Eq. (jlj) we infer that the link 
operators Kj, Kg, and act on u as X^, Y^, and up to a permutation. Conjugating 
Kf, Kg, and K^ by an appropriate single-qubit Clifford unitary operator Vu one can 
generate any permutation of Yu, and In addition, the depolarizing noise model 
used in this paper assigns equal probabilities to all Pauli errors on u, and thus it is 
also invariant under permutations of Xu, Yu, and Zu- It shows that the specific form of 
the link operators does not play any role as long as they obey the commutation rules 
Eq. (jlj). Let Q = {K^, e G E) he the gauge group generated by the link operators. 

Definition 1. A subset of links M (1 E is called a cycle iff every site has an even 
number of incident links from M. 

Lemma 3.1. The subsystem code Q has no logical qubits, that is, C{Q) C Q. 
Furthermore, P G C{Q) iff P ^ IleeM -^e for some cycle M 'O E. 

Proof. It suffices to prove the second statement of the lemma. Consider any operator 
P G C{Q). We shall say that P locally anticommutes with a link operator Ku^v iff 

P T<^ — — K P nnH P T<^ — —K P 

The subset M O E including all links locally anticommuting with P G C{Q) will be 
called an anticommuting set of P. We claim that the anticommuting set M is a cycle, 
that is, any site has even number of incident links from M. Indeed, Eq. (jlj) implies that 
for any site u the product of K^. over the three links e incident to u acts trivially on u. 
Hence P can locally anticommute only with even number of Ke among the three links 
e incident to u. It shows that M is a cycle. Consider now an operator 

eGM 

Using the fact that M is a cycle and Eq. (jlj), one can easily check that P' G C{Q). 
Furthermore, the anticommuting set of P' is M. Hence PP' has trivial anticommuting 
set, that is, PP' locally commutes with any link operator Kg. It means that PP' 
commutes with any single-qubit Pauli operator, which is possible only if PP' ~ /. 
Hence P ~ P'. □ 



3.3. Loop Operators on 3-valent Hypergraphs 

As we demonstrated in the previous section, the honeycomb model and its 
generalizations to 3-valent graphs give rise to trivial subsystem codes encoding no logical 
qubits. In particular, any pair of loop operators commute with each other regardless 
of the homological class of the loops. Hence the codes derived from the honeycomb 
model cannot satisfy condition (C3). To overcome this problem, it is helpful to change 
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the perspective and view as fundamental objects the loop operators rather than the 
link operators. Using the ideas from the previous section, we shall define a group of 
loop operators Qioop associated with any 3-valent hypergraph £ = (V, E) in which edges 
represent pairs and triples of sites. The latter type of edges gives rise to anticommuting 
pairs of loop operators which could be used to encode logical qubits. We shall see 
that a pair of loop operators anticommute (commute) iff the two loops share an odd 
(even) number of triangles, see Lemma [3?2l Once the group of loop operators is defined, 
the gauge group ^ of a subsystem code and its stabilizer group S = Q \~\ C{Q) can be 
reconstructed using identities 



We will prove that Q has two-qubit generators for all hypergraphs C obeying certain mild 
restrictions. The main advantage of this top-bottom approach is that once a suitable 
group of loop operators is chosen, we can easily find bare logical operators Xj, Zi and 
generators of the stabilizer group S. Indeed, Eq. ([71) implies that 



that is, the bare logical operators Xi,Zi G C{Q) can be chosen as pairs of loop 
operators obeying the standard Pauli commutation rules, while the stabilizers are the 
loop operators commuting with all other loop operators. In contrast, deriving logical 
operators and stabilizers starting from the link operators may require rather lengthy 
calculations and does not give much insight into the topological nature of the code. 

As before, we consider a system of n physical qubits occupying a set of sites V 
embedded into a plane or a torus. To define a subsystem code, we shall identify V with 
a set of vertices of a hypergraph C = (V, E) with two species of edges: edges connecting 
pairs of sites and edges connecting triples of sites. For consistency with the previous 
notations, we shall refer to the edges of the first type as links and the edges of the 
second type as triangles. Thus E = E2U E^, where E2 is a set of links and E^ is a 
set of triangles. More specifically, a link {u, v) G E2 is an unordered pair of sites and 
a triangle {u,v,w) G E3 is an unordered triple of sites. We shall say that an edge e is 
incident to a site u iff e = {u,v) or e = {u,v,w). There will be only three restrictions 
imposed on the hypergraph £, namely, 

• Each site has exactly three incident edges, 

• Any pair of edges share at most one site, 

• All triangles are pairwise disjoint. 

We show examples of such hypergraphs C in Figures [2l[5l A different example is the 
triangle-honeycomb lattice studied in Ref. pB] but for this decorated honeycomb model 
Gioop is trivial. By analogy with the honeycomb model, we consider a family of edge 
operators Kf. living on edges e G -E and satisfying the commutation rules Eq. (jl]). We 
shall refer to operators as link operators if e G £"2 and as triangle operators if e G £^3. 
The same arguments as above show that the edge operators are uniquely defined by the 
commutation rules Eq. (jlj) up to local unitaries. Also the commutation rules Eq. (jlj) and 



Q = C{Qioop) and S = Qioop H C{Qioop)- 



(7) 



C{G) = Qioop, 



(8) 
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our assumptions on C imply that the triangle operators are three-qubit Pauli operators, 
for example, K^vw = Z^Z^Z^. 



Definition 2. A subset of edges M C is called a cycle iff every site has even number 
of incident edges from M. 

Our assumptions on C imply that any cycle M consist of a disjoint union of loops, 
open paths, and triangles such that free ends of every open path belong to some triangle, 
and every triangle {u,v,w) G M has exactly three open paths terminating at u,v, and 
w. See Figures I2]|5] for examples of a cycle. Given any cycle M, the operator 

W{M) =11 (9) 

eeM 

will be referred to as a loop operator associated with M. In contrast to the standard 
honeycomb model of Section [321 the loop operators W{M) may anticommute with each 
other. 

Lemma 3.2. Let e & E be any edge and M (1 E be any cycle. The loop operator W {M) 
anticommutes with Ke if e is a triangle contained in M , i.e, e G Mn£^3, and commutes 
otherwise. Furthermore, for any cycles M, M' C E one has 

W{M)W{M') = (-i)^W^') W{M')W{M), (10) 

where A(M, M') is the number of triangles shared by M and M' . 

Proof. Let Au{e) be the incidence matrix of £, that is, Au{e) = 1 if e is an edge incident 
to u and Au{e) = otherwise. Recall that K^Ke' = (— 1)'''('^''^') K^'K^, where //(e, e') = 1 
iff e and e' share exactly one site and T]{e, e') = otherwise. One can easily check that 
for any pair of edges e, e' G -E one has the following identity: 

vie, e') = X3(e)4,e' + E Me') (mod 2), (11) 

where Xsi^) = 1 if e is a triangle and Xsi^) = if e is a link. For any site u and for any 
cycle M one has 

E^.(e')=0 (mod 2). 

e'gAf 

Taking into account Eq. ffTTl) one arrives at K^W{M) = (-l)°(^)iy(M)ire, where 
(^(e) = E v{e,e') (mod 2) 

e'eAf 

= E X3(e)5e,e' + E^"(e') (mod 2) 

e'GAf uee 

= X3(e) E ^e.e' ("^Od 2). 
e'eAf 

Hence K^. anticommutes with W{M) if e G M is a triangle and commutes otherwise. It 
immediately implies Eq. ( ITOj) . □ 



Constructions and noise threshold of topological subsystem codes 



11 



It is worth pointing out that the number of triangles in any cycle must be even, 
so Eq. ffTOj) is consistent with the fact that W{M) commutes with itself. Indeed, each 
open path in M contributes two sites that must be covered by triangles in M. Since 
each triangle covers exactly three sites, the number of triangles must be even. 

3.4- Reconstructing the Gauge Group from the Loop Operators 

Recall that our goal is to represent the group of loop operators Qioop as a centralizer of 
the gauge group Q which has two-qubit generators, see Eqs. (ITfSj) . Let us first explain 
how to transform the hypergraph C = (V, E) into an ordinary graph C = (V, E') which 
will be the underlying lattice for our subsystem code. The graph C has the same set of 
sites V representing locations of the physical qubits. Recall that E = E2U E^, where 
E2 is the set of links and E^ is the set of triangles. Every link [u, v) G E2 becomes a 
link of the graph C, that is, we have the inclusion E2 ^ E' . Furthermore, every triangle 
(m, V, w) G -E3 contributes a triple of links to the graph namely the links (m, f ), {y, w), 
and {u, w). In other words, the graph C is obtained from the hypergraph C by breaking 
up all triangles into the individual links. To avoid confusions between various species 
of links and edges, we shall refer to the links of C originating from links of C as solid 
links. The links of £' originating from triangles of C will be referred to as dashed links. 
Hence E' = EsU Ed where Eg = E2 and Ed are the sets of solid and dashed links. The 
solid and dashed links will be rendered properly on all figures. 

Now we are ready to define the gauge group Q. As was mentioned above, the 
commutation rules Eq. (j4]) uniquely specify the edge operators K^. up to local unitaries. 
It will be convenient to fix the choice of these local unitaries such that all triangle 
operators include only Pauli Z, that is, 

Ku,v,w = ZuZ^Z^ for all G £"3. (12) 

We define the set of link operators K'^, e E E' generating the gauge group Q as follows: 

„ = Ku,v if (m, f) is a solid link, 
K'^ J, = ZuZ^ if (m, v) is a dashed link. 

In other words, each triangle operator Ku.v.w gives rise to a triple of dashed link operators 
obtained by restricting Ku^^^w on pairs of qubits iy^w), and {u,w). (Note that 

these dashed link operators are not independent since K'^^K'^^K'^^ = I). The main 
result of this section is the following lemma. 

Lemma 3.3. Let Q be the group generated by solid and dashed link operators K'^, e E E' . 
Then Qioop = C{Q). 

Proof. Let us first prove the inclusion Qioop ^ C{Q). We need to check that for any link 
e G -E' and for any cycle M in the hypergraph C the link operator K'^ commutes with 
the loop operator W{M). 

Suppose e G -Es is a solid link. Then K'^ = Kg. Since e is not a triangle, Ke 
commutes with W{M), see Lemma [3.21 
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Suppose e = {u,v) G is a dashed link so that K^^^ = Z^Z^. Let / = 
(u, V, w) G -E3 be the corresponding triangle in C. Using the assumption Eq. f|T2l) 
and the commutation rules Eq. (jl]), we conclude that each of the sites u,v,w has two 
incident links in £ and the corresponding link operators act on these sites by X and Y. 
li f ^ M, each site u, v, w has either zero or two incident links from M. Hence W{M) 
acts on u,v,w by / or Z. Since K'^ = Z^Z^, it follows that K'^ commutes with W{M). 
If / G M, each site u, v, w has exactly one incident link from M. Hence W{M) acts on 
u,v,w by X or Y. Again we conclude that K'^ = Z^Z^ commutes with W{M). 

It remains to prove the reverse inclusion C(^) C Qi^oji- The proof of this part is 
almost identical to the proof of Lemma ISTTl Consider any operator P G C(^). We claim 
that for any edge operator i^e, one has one of the two possibilities: (i) P commutes 
with K(. at every qubit in the support of i^e, or (ii) P anticommutes with K^. at every 
qubit in the support of K^. Indeed, if e = [u^v] G -^2 is a link then Ke = K'^ E Q 
is a two-qubit Pauli operator commuting with P, that is, we have either (i) or (ii). If 
e = (m, V, w) G -E3 is a triangle, the assumption P G C{Q) means that P commutes with 
ZuZy, Z^Zyj, and ZuZ.^^. Since = Z^Z^Z^, we have either (i) or (ii). We shall say 
that P locally commutes with Kg in the case (i) and locally anticommutes with Ke in the 
case (ii). The subset M ^ E including all edges locally anticommuting with P G C{Q) 
will be called the anticommuting set of P. 

We claim that the anticommuting set M is a cycle, that is, any site u has even 
number of incident edges from M. Indeed, consider any site u and let f, g, h E E he the 
triple of edges incident to u. Using Eq. (jl]) and our assumptions on the hypergraph C, 
we conclude that Kf, Kg, acts on u as X, F, Z up to a permutation. It follows that 
Pu anticommutes with even number of operators Kj, Kg, Kh- The observation made 
above shows that P locally anticommutes with even number of Kj, Kg, K^, that is, M 
is a cycle. 

Consider a loop operator W{M) = YieeM Ke G Gioop- The inclusion Qioop ^ C(Q) 
implies W{M) G C(^). A direct inspection shows that the anticommuting set of W{M) 
is M. Hence P-W{M) has trivial anticommuting set, that is, P-W{M) locally commutes 
with any edge operator K^. It means that P ■ W{M) commutes with any single-qubit 
Pauh operator, that is, P ■ W{M) ~ J and P ~ W{M). Hence P G Gioop- □ 

3.5. Examples of 3-valent Hypergraphs 

In this section, we consider several explicit examples of 3-valent hypergraphs and find 
the corresponding group of loop operators Qioop- As was argued above, the group Qioop 
alone determines all properties of the code, so we shall not need the explicit form of the 
gauge generators. 

Throughout this section we shall use a notation Z{C) for the binary cycle space of 
a hypergraph C = {V, E). In other words, elements of Z{C) are cycles M O E and the 
addition of cycles corresponds to the symmetric difference. Clearly, choosing generators 
of the group Qioop is equivalent to choosing basis vectors of the binary cycle space Z{C). 
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Let us first consider the hypergraph Cgo corresponding to the square-octagon lattice 
in Figure [2l Every site of Cgo has two incident hnks (solid hues) and one incident 
triangle (filled triangles). Let us first list basis cycles that have spatially local support 
(homologically trivial loops). Such cycles are associated with the elementary octagons 
and squares. Each octagon contributes two basis cycles A, C e Z{Cso), see Figure O 
The cycle A consists of 8 links forming the boundary of the octagon. The cycle C is a 
loop of 8 triangles and 12 links encircling the octagon (the links and triangles forming 
each cycle are highlighted in red in Figure El) . Analogously, each square contributes two 
basis cycles B,D E Z{Cso)- The cycle D consists of four links forming the boundary of 
the square. The cycle 5 is a loop of 4 triangles and 6 links encircling the square. Let us 
consider now the corresponding loop operators, for example, W{C) = Yie&c ^e- Note 
that each site u has either or 2 incident edges belonging to C. Commutation rules 
Eq. (jll) imply that if u has two incident edges from C, the corresponding edge operators 
Kf. act on u by different Pauli operators. It means that W{C) acts non-trivially on a 
qubit u iS u has two incident edges from C. A direct inspection shows that W{C) acts 
non-trivially on 24 qubits. Similar arguments show that W{A),W{B), and W{D) act 
non-trivially on 8, 12, and 4 qubits. 

The remaining basis cycles are those associated with the homologically non-trivial 
loops on the lattice. We denote these cycles Zi, Z2,Xi, and X2 see Figure [3] and Figure HI 
A simple algebra shows that the cycles A, B, C, D and Zi, Z2, Xi, X2 form a basis of the 
cycle space Z{Cso)- 

Let us now consider commutation rules between the corresponding loop operators. 
We note that the cycles A and D contain no triangles which implies that 14^ (A) and 




Figure 2. (Color Online) The 3-valent hypergraph Cso — {V,E) associated with 
the square-octagon lattice embedded into a torus (periodic boundary conditions) . The 
hypergraph has two types of edges: links (solid lines) and triangles (gray-filled and red- 
filled triangles) . Each site has two incident links and one incident triangle. A cycle is a 
subset of edges MCE such that each site has even number of incident edges from M. 
The hypergraph has four types of elementary cycles indicated by letters A, B, C, D. 
The links and triangles forming each cycle are highlighted in red. The underlying 
lattice of the subsystem code is obtained from the hypergraph Cso by breaking up 
each triangle into a triple of links (dashed lines) . 
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Figure 3. (Color Online) Basis cycles Zi, Z2 G Z{Cso) corresponding to homologically 
non-trivial vertical loops on the lattice. 




Figure 4. (Color Online) Basis cycles Xi,X2 E Z{Cso) corresponding to 
homologically non-trivial horizontal loops on the lattice. 

W{D) commute with all loop operators, see Lemma 13. 2[ that is, W{A),W{D) e 
Gioop n C{Qioop)- Since Qioop n CiGioop) = S, see Eq. (j?!), we conclude that W{A) 
and W{D) are stabilizers. Furthermore, one can easily check that any cycle of type 
B OT C shares either zero or two triangles with any other cycle B,C, Zi, Z2,Xi,X2. 
Lemma [3.21 then implies that W{B) and W{C) commute with all loop operators, that 
is, W{B), W{C) E S. Thus all the loop operators associated with A, B, C, and D cycles 
are stabilizers. 

The cycles Zi, Z2 share exactly one triangle with the cycles Xi, X2 respectively, see 
Figures |31H1 Lemma 13.21 implies that 

W{X,)W{Z,) = {-lY-'^W{Z,)W{Xa), a, 6 = 1,2. (14) 

We conclude that the code Q = C{Qioop) has exactly two logical qubits and the loop 
operators W{Zi),W{Z2), W{Xi),W{X2) can be chosen as the bare logical operators 
Zi, Z2, Xi, X2 respectively. 

Let us mention that the subsystem code Q = C{Qioop) coincides with the topological 
subsystem code of [9] associated with the Union Jack lattice (the dual of the square- 
octagon lattice). Hence applying our construction to the square-octagon lattice does not 
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provide new examples of codes but rather provides a new angle to look at the already 
known code. 

Let us now consider the hypergraph C^s corresponding to the lattice shown in 
Figure [5l We shall refer to this lattice as the five-squares lattice since it unit cell 
consists of five squares. In contrast to the square-octagon lattice, some sites of £5^ 
have three incident links and others have two incident links and one incident triangle. 
Simple algebra shows that the cycle space Z^C^g) has four types of local basis cycles 
A, B, C, D associated with loops encircling squares and octagons, see Figure [5l and 
homologically non-trivial basis cycles Zi, Z2, Xi, X2, see Figure |6] and Figure [71 As 
before. Lemma [32] implies that the loop operators W{B), W{C), W{D) generate 

the stabilizer group S = Gioop H C{Qioop)- The loop operators associated with cycles 
Zi, Z2,Xi,X2 obey commutation rules Eq. ( !T4l) and thus can be chosen as bare logical 
operators on the two logical qubits. 

Let us briefly sketch how the syndrome measurement for the stabilizers of different 
type can be simulated by measuring eigenvalue of the two-qubit link operators K'^. 
Stabilizers of type A and D correspond to closed even-length loops containing no 
triangles. Hence the syndrome measurement strategy used for the honeycomb model 
in Section [3l2] can be applied directly to stabilizers of type A and D. Consider now a 
stabilizer of type B. It consists of a solid center square and a length-8 external loop 
with alternating solid and dashed links, see Figures [21151 One can easily check that the 
stabilizer W{B) can be represented as a product of the eight link operators K'^ lying on 
the external loop and operators K'^ associated with the two vertical solid links e ^ B 
lying on the central square, see Figure [HI In order to satisfy condition Eq. one can 
divide these 10 links into three levels: Level-1: horizontal solid links, Level-2: vertical 
solid links, and Level-3: dashed links. One can easily check that all link operators K'^ 
at the same level mutually commute and any link operator at level j commutes with 




Figure 5. (Color Online) The 3-valent hypergraph = associated with 

the five-squares lattice embedded into a torus (periodic boundary conditions). The 
underlying lattice of the subsystem code is obtained from the hypergraph £5^ by 
breaking up each triangle into a triple of links (dashed lines). Four type of local 
basis cycles A, B,C,D £ Z{Cso) are highlighted in red. 
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Figure 6. (Color Online) Basis cycles Zi, Z2 £ Z^C^s) corresponding to homologically 
non-trivial vertical loops on the lattice. 

,R A ,n A ,n A 

anitf Vi;:p:r/\j:;:n:r/ 

Figure 7. (Color Online) Basis cycles Xi,X2 € ^(£5^) corresponding to 
homologically non-trivial horizontal loops on the lattice. 

the product of all link operators at levels < j. Hence we can satisfy Eq. ([6]) using any 
"level-nondecreasing" order of measurements: we first measure all level- 1 links (in any 
order), then measure all level-2 links (in any order), and finally measure all level-3 links 
(in any order). A similar strategy can be applied to stabilizers of type C. For example, 
the stabilizer W{C) for the five-squares code, see Figure [5l can be represented as a 
product of 40 links operators K'^. These 40 links can be divided into three levels such 
that all link operators K'^ at the same level mutually commute and any link operator 
at level j commutes with the product of all link operators at levels < j. Hence any 
level-nondecreasing order of measurements obeys the condition Eq. ([6]). 

4. Decoding Algorithms 

4-1. Maximum Likelihood Decoding 

Here we describe the maximum likelihood decoding for an arbitrary subsystem code Q. 
Our derivation mostly follows [9]. Let us assume that the noise can be described by 
random Pauli errors E drawn from some distribution Pr [E). Let S = Q r\C{Q) he the 




Constructions and noise threshold of topological subsystem codes 



17 



stabilizer group of Q. For any Pauli error E let [E] be the (left) coset of Q that contains 
E, that is, [E] = E ■ Q. Obviously, the syndrome s caused by an error E depends only 
on the coset [E]. Let s{h) be the syndrome caused by any error in the coset h. The goal 
of the decoding algorithm is to identify an error up to a gauge operator based on its 
syndrome. In other words, one needs to guess the coset h based on the syndrome s{h). 

Proposition 1. Suppose the code has k logical qubits. Then for any given syndrome s 
there are 4*^ different cosets h such that s{h) = s. 

Proof. Indeed, suppose h = [E] and h' = [E'] such that s{h) = s{h') = s. Then 
E-E' e C{S). Using the identity C{S) = C{g) -Qwc conclude that E-E' = PG for some 
bare logical operator P G C{Q) (may be the identity), and some gauge operator G E Q. 
Since E and E' arc defined modulo gauge operators, we can ignore G. Furthermore, 
if h and h' are distinct cosets then P is a non-trivial logical operator. Since there are 
exactly 4*^ distinct bare logical operators, there must be exactly 4*^ distinct cosets h such 
that s{h) — s. □ 

Let S{s) be the set of all 4^ cosets h such that s{h) = s. The chosen probability 
distribution of errors induces some probability distribution P{h) of cosets, namely, 

p{h) = E (E) ■ 

Eeh 

Given a syndrome s, the maximum likelihood decoder computes the probability P{h) 
for all cosets h e £{s) and then selects the coset with the largest probability. For any 
error E one has 

Pr([^])=EP^(^G)- 
Gee 

For example, for depolarizing noise with an error rate p, one has 

/ \ \m 

Pr{E) ' ^ > 



,3(1 -p)y 

Let us introduce an inverse temperature (3 such that 

Then one arrives at 

Pr([^])- Y.e-^^\^''\ = Z{E). 

Note that by definition the partition function Z{E) depends only on the coset [E]. Here 
Z{E) can be regarded as a partition function of a classical spin Hamiltonian with a 
quenched disorder determined by E and the inverse temperature 2^. The classical spins 
are associated with generators K^. of the gauge group. Indeed, we can represent 

G^l[Kf^\ a(e)e{0,l}. 

Then the Hamming weight \EG\ can be expressed a local Hamiltonian depending on the 
classical binary variables a(e). Parameterizing a's by Ising spins, one express Z{E) as 
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a partition function of an Ising-like Hamiltonian with inverse temperature /3. Finding 
the values of the binary variables a{e) which minimize the weight \EG\ corresponds to 
minimizing the energy of the local Hamiltonian. Due to the particular structure of the 
Pauli group, this minimization is an instance of a MAX-XORSAT problem with local 
SAT constraints between binary variables on a 2D lattice. 

Following [71 [9] we can relate the threshold value of the error rate p to a phase 
transition in the statistical mechanics model corresponding to Z{E), where E is regarded 
as a quenched disorder. Indeed, suppose the error rate p is below the threshold and let 
E he a typical error with a syndrome s = s{E). Then one should expect that in the 
limit of large lattice size L all cosets in S{s) have negligible probability except for the 
coset [E]. In other words, Pr([i?]) ^ Ft (jPE]j for any non-trivial logical operator 
P € C{Q)\Q. On the other hand, if the error rate p is above the threshold, one should 
expect that in the limit of large lattice size all cosets in S{s) have the same probability. 
Thus we can determine the threshold error rate Pc using the following heuristic criterion: 

p < Pc =^ lim Pr (E) log 

P>P. ^ limEPr(i?)log(|H^)=0. (16) 

Assuming that the code Q has a macroscopic distance, the logical operator P can be 
"cleaned out" from any local region of the lattice by a gauge transformation [29]. Hence 
the difference between the partition functions Z{E) and Z{PE) is hidden if one inspects 
any local region of the lattice. It can only reveal itself if one inspects the entire lattice 
and if the system exhibits a long-range order. Hence one should expect that the point 
p = Pc corresponds to a phase transition between a magnetically ordered phase {p < Pc) 
and a paramagnetic phase {p > pc). 

In general, implementing maximum likelihood decoding is a computationally 
inefficient task, since one would have to estimate the different partition functions Z{PE) 
for a local 2D Hamiltonian with quenched disorder. Thus one has to seek efficient 
strategies which approximate the action of such ideal decoder. Given the relatively low 
error rates p which correspond to low temperatures via Eq. (ITSll . one may argue (as 
in [7]) that maximum-likelihood decoding is pretty well approximated by considering 
/3 — ?► oo or doing minimum-weight decoding. For minimum-weight decoding, we consider 
which 2D local Hamiltonian with quenched disorder determined by PE has minimum 
energy. In other words, we find the error consistent with the syndrome which has 
minimum Hamming weight. Again, for this minimization problem, -instances of a two- 
dimensional MAX-XORSAT problem-, there does not necessarily exist an exact efficient 
algorithm. For the 2D surface code, there is an efficient algorithm for minimum-weight 
decoding which is Edmond's efficient polynomial time (in the lattice size) minimum- 
weight matching algorithm [30] . Of course, even the absence of an exact algorithm does 
not preclude the existence of good decoders based on heuristic optimization strategies 
which may be more computationally efficient than exact algorithms. This is the 
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motivation for the work in [TU| [ST] in which the authors develop fast renormahzation 
group decoding algorithms which are suitable for any topological stabilizer or subsystem 
code. 

In Sections 14.21 and 14.31 we will describe a simple decoding algorithm for the five- 
squares code. The decoding algorithm in Section 14.21 is improved in Section 14.31 by 
reducing the weight of the guessed error which gives a slight boost to the noise threshold. 
For the five-squares code, we have developed a decoder which is fine-tuned to the 
particular structure of the code and which maps the decoding after some preprocessing 
onto decoding of two copies of the toric code. For the decoding of the toric code, one 
can then use a minimum-weight matching algorithm or renormahzation group decoding. 
Let us point out that applying a similar decoding strategy to the subsystem code defined 
on the square-octagon lattice, see Figure O one can reduce the decoding to correcting 
Z errors using stabilizers of the topological color code [32]. In |3l] the authors are 
able to establish a much more general correspondence between the toric code and any 
topological subsystem code [9]. They show that any topological subsystem code is locally 
unitarily equivalent to some copies of the toric code. This implies that the decoding 
of any topological subsystem code can proceed via efficient decoding algorithms for the 
toric code. 

4-2. Decoding for the Five- Squares Code 

To describe our decoding algorithm it will be convenient to fix a particular choice of the 
link operators associated with the solid links. We will choose Ku^v = XyY.f, for all 
solid links (m, v) that belong to the boundary of some square, see FigurelHl Commutation 
rules Eq. (j4]) then imply that one has to choose Ku,^ = Z^Z^ whenever {u,v) is a solid 
link such that u and v belong to different squares, see Figure [91 Recall that the link 
operators associated with dashed links are fixed by Eq. ( 1131) . 

Our decoding algorithm exploits the fact that errors only need to be corrected 
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Figure 8. Link operators Ku,v associated with, solid links {u, v) that belong to 
boundary of some square. 



Figure 9. Link operators Ku^v associated with solid links {u, v) connecting two 
different squares. 



Constructions and noise threshold of topological subsystem codes 



20 



2 

4i [3 

/ I I ^ 
5 6 19 101 ^sl3 14 



/ 1 I 
^ 19 Id 

n;n:D 

8 7 N 112 111 ,^16 i; 

□ ::n 



(a) It is sufBcient to (b) Z errors need to be 
correct X errors at qubit corrected at qubits 1,2, 
1. 4, 19, and 20. 

Figure 10. Locations of the errors. 



up to multiplication with elements in the gauge group Q. In particular, gauge group 
elements can be used to move errors around and simplify the decoding. To illustrate this 
point let us consider any solid square with the corresponding qubits labeled by 1, 2, 3, 4, 



see Figure 10(a) We claim that instead of considering general Pauli errors acting on 
these qubits, it suffices to consider simplified errors from the group (Xi, Zi, Z2, Z3, Z4), 
that is, arbitrary Z errors and X errors acting only on qubit 1. Indeed, using the 
gauge operators Ku,v = XuY^ and Z errors one can generate any Pauli operator P that 
contains Z's and even number of X's. Multiplying P by Xi we can obtain any Pauli 
operator on qubits 1, 2, 3, 4. Thus it suffices to consider arbitrary Z errors and X errors 
acting only at a single location in each square (say, the north-west corner). Let us also 
note that the D-stabilizer associated with the considered square is 

We can now use a similar trick to reduce the set of possible locations for Z errors. 
Consider a unit cell of the lattice shown in Figure 10(b) that consists of 20 qubits. We 
claim that any Z error acting on the unit cell can be reduced to Z errors at the locations 
1,2,4,19,20. Indeed, let us first describe relevant gauge operators of Z type that we 
shall use to simplify the error. Define a group C (Zi, . . . , Z20) generated by the link 
operators ZuZ^ associated with 12 dashed links within the unit cell and D-stabilizers 
ZsZtZuZ^ associated with the four squares on the top, bottom, left, and right edges 
of the cell. Let us partition the unit cell into 12 internal qubits that belong to some 



triangle and 8 external qubits, that is, qubits 1,2,5,8,14,15,19,20, see Figure 10(b) 
Consider a subgroup Qf^^ generated by the dashed link operators ZuZ^, and truncated 
versions of the D-stabilizers .^3^4, Z^Z^, Zi^Ziq, and ZiyZig. In other words, is a 
'projection' of the group onto the internal qubits. Since any pair of internal qubits 
can be connected by a path of ZZ generators belonging to Q-^^, we conclude that ^f^^ 
contains all even-weight Z-operators on the internal qubits. It follows that any Z error 
P acting on the unit cell can be represented as P = P'G, where G & and P' contains 
only Z errors at the external qubits and some selected internal qubit, say qubit 4. 
Finally, we can move Z errors out of any external qubit u G {5, 8, 14, 15} using the solid 
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link operators Z^Z^, where v G {1, 2, 19, 20} belongs to one of the adjacent unit cells. 

To summarize, it suffices to consider only X errors acting on the north-west corners 
of each square and Z errors acting on locations 1,2,4,19,20 of each unit cell, see 
Figures [IO(a)]flO(b) 



We choose the following decoding procedure to correct this reduced set of errors. 
Assume that all stabilizer generators A, B and C and D have been measured perfectly. 
We first consider the syndromes of stabilizers D. Note that D anti-commutes with any 
single X error at the corresponding square and commutes with all Z errors. For all 
non-zero D-syndromes, we will apply X to the north-west corner of the corresponding 
square which corrects all X errors. After this step we are left with Z errors at the 
locations {1, 2, 4, 19, 20} in each unit cell. 

Note that this error correction step may change the syndromes of A, B, and C. 
For example, an X error in the square in the center of the unit cell anticommutes with 
the corresponding stabilizer B. The syndromes A, B, and C can be either measured 
after the first error correction step, or their new values can be calculated based on their 
initial values and the error corrections that were applied (i.e. updating the so-called 
Pauli frame). 

Any B stabilizer acts trivially on locations {1,2, 19,20} in the corresponding unit 
cell and acts by Y at the location 4. Hence the only Z error in the unit cell that can 
create a non-trivial i?-syndrome is the error Z4. For all non-zero i?-syndromes we will 
apply Z4 in the corresponding unit cell to correct the Z-error. 

It remains to show how to correct the remaining Z errors at locations {1, 2, 19, 20}. 
These qubits are divided into two groups, which are depicted separately in Figure [11] 
and [121 These figures show how the stabilizers A and C act on the qubits. Each 
of the stabilizers acts on four of the considered qubits. Stabilizers of type A are 
represented by blue and orange squares, and C operators are represented by red and 
green squares. Thus the reduced lattice effectively consists of two disjoint lattices on 
which the stabilizers acts as four-body plaquette operators. The only occurring errors 
on these two disjoint lattice are Z errors which anti-commute individually with the A 
and C stabilizers which cover their support. Hence, one can view the lattice in e.g. 
Figure [H] as half of a toric code and the Z errors can be corrected using any decoding 
strategy for the toric code. We have used decoding via the minimum weight perfect 
matching algorithm similarly as in [7]. 

We can now try to give a rough estimate of the noiseless error correction threshold, 
solely based on this mapping onto decoding for the toric code. An upper-bound for the 
threshold against, say, X errors for the toric code is 11% [20] and this is close that what 
minimum weight decoding achieves. We can derive a rough estimate for the effective 
single-qubit noise rate for Z errors on the reduced toric code lattice, as a function of the 
basic depolarizing noise rate. On the original lattice, we have a single-qubit Z-error rate 
of 2p/3 where p is the depolarizing error rate. Some X errors get corrected, but some 
also get mapped onto Z errors so, perhaps one should consider that after the first round 
of error correction (the processing of the syndrome of operator D), one has Z error rates 
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Figure 11. (Color Online) The stabilizers A and C on a lattice with the four Z errors 
in the unit cell. 
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Figure 12. (Color Online) The stabilizers A and C on a lattice with the four Z errors 
in the unit cell. 



between 2p/3 and p on all sites of the lattices (for example, the Z error rate on qubit at 
the location "1" is ^p/Vj. Using multiplication by the gauge group elements, we replace 
these 20 possible error locations in a unit cell by 4 effective error locations. Hence the 
effective Z-error on these 4 error locations is between lOp/3 and bp. This would suggest 
that the five-squares code can achieve at most a depolarizing noise threshold between 
2.2% and 3.3% using this decoding method. 

This analysis is rough, but it suggests why this decoding strategy for the fives- 



Constructions and noise threshold of topological subsystem codes 



23 



squares code will lead to a threshold which is markedly lower than the threshold of the 
surface code. The topological subsystem code has a larger overhead in terms of number 
of qubits in order to make for very local gauge group generators. If this larger overhead 
does not contribute to the error correction, it simply increases the effective error rate 
that is seen by the toric code decoder. 

In this rough estimate we have also neglected correlations between errors. As we 
mentioned earlier, elementary single qubit errors can be mapped onto two-qubit errors 
by multiplication with gauge group elements. Hence the effective error model that is 
seen by the toric code decoder is not an i.i.d. single-qubit error model. Our experimental 
determination of the threshold simulation of course takes these correlations into account. 
We believe that correlations do not play a major role in depressing the threshold. 



4-3. Improved Decoding for the Five-Squares Code 

The decoding algorithm in the previous subsection is suboptimal because it does not 
minimize the weight of the corrected error. We reduce the weight of the corrected error 
by modifying our algorithm - we correct the X errors in the first step of the algorithm 
at such locations so that we reduce the number of Z errors that result. 

First, we observe that the X errors do not necessarily need to be corrected at 



qubit 1 as indicated in Figure 10(a) Instead, if the D syndrome corresponding to a 
particular square is nontrivial, we can correct the concerned bit flip at any single corner 
of the square. Further, we observe that depending on the location where we perform 
the correction, the B syndrome in the corresponding unit cell will or will not change its 
sign. This is illustrated in Figure [I3l 

To describe the main correction step, let's focus on a particular unit cell and its 
syndromes B and D. Let's first assume that at least one of the five D syndromes is 
non-trivial, and the B syndrome is trivial. Then, if we correct the bit flips indicated by 
the D syndromes at the locations that do not affect the B syndrome, we know that the 
Z error will not occur at qubit 4. If there is a least one non-trivial D syndrome, and the 
B syndrome is non-trivial as well, we perform exactly one of the X error corrections at 
a location that changes the B syndrome. Again, the Z error will not occur at qubit 4. 
The last case to analyze has all D syndromes trivial, and the B syndrome non-trivial. 
In such a case the Z error correction at qubit 4 still needs to be performed, and the 
improved algorithm does not reduce the weight of the corrected error. In the first two 
analyzed cases, the weight of the corrected error decreases by up to one. 



4-4- Syndrome Measurements in the Five-Squares Code 

An important aspect concerning the practicality of the five-squares code, in particular 
the noise threshold against noisy error correction, concerns the number of two-qubit 
gauge element measurement needed to determine the syndrome of the stabilizer 
generators. The probability with which a stabilizer generator has a faulty syndrome 
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Figure 13. Illustration of the improved decoding. 

Figure 14. Example - these ten link operators need to be measured to determine 
syndrome B of the five-squares code. 



is experimentally determined by the error rate on the two-qubit gauge element 
measurements. 

Consider, for example, stabilizer B in the five-squares code. Figure [H] shows that 
10 link operators need to be measured to determine the syndrome. It is easy to verify 
that to determine syndromes A, C, and D in the five-squares code, respectively, 8, 
40, and 4 two-qubit measurements need to be taken. One may then anticipate that 
the probability of a faulty stabilizer generator syndromes, say for a stabilizer of type 
C is roughly at most 40 times the probability for a faulty two-qubit gauge element 
measurement. This implies that the noise threshold against noisy error correction may 
be considerably lower than the 2% noise rate that we obtain here. We note that for 
the square-octagon code the number of measurements to determine the four types of 
syndromes is slightly smaller, namely 14, 8, 4, and 20. 

5. Experimental Evaluation of the Threshold 

We performed extensive simulations to evaluate the threshold properties of the simple 
and improved decoding algorithm. We experimented with a range of lattice sizes and 
depolarizing error rates, recording the ratio of unsuccessful error corrections for each 
such configuration. This allowed us to accurately estimate the threshold because for 
large lattice sizes, the decoding succeeds with high probability for error rates below the 
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threshold, and fails otherwise. 

Our simulation first generates a random depolarizing error E, and then it utilizes 
the simple or improved decoding algorithm to determine the error Eguessed that should be 
corrected according to the decoder. In the final step, we determine whether correcting 
the original error E by applying the guessed error Eguessed fails, i.e., if EEguessed is 
proportional to a logical operator. More precisely, for a particular depolarizing error 
rate p, and a fixed lattice size, the simulation proceeds as follows: 

(i) Generate a random error E with depolarizing probability p in the five-squares 
lattice. 

(ii) Calculate syndromes A, B, C, and D. 

(iii) Add X errors indicated by syndromes D to Eguessed, and update syndromes A, 
B, and C assuming that the error Eguessed is corrected (i.e., the Pauli frame is 
updated). 

(iv) Add Z errors indicated by syndromes B to Eguessed, and update syndromes A, and 
C accordingly. 

(v) Correct the errors in the two sub-lattices in Figure [H] and [12] separately. Use 
the minimum weight matching algorithm, and add Z errors on the shortest lines 
connecting the matched non-trivial syndromes to Eguessed, as described in detail 
below. 

(vi) Calculate O = EEguessed- If O is a logical operator on the five-squares lattice, 
correction fails. Otherwise it succeeds. To determine if O is a logical operator, we 
only need to check if it anticommutes with one of the logical operators Xi or Zi. 

(vii) The simulation above is repeated for a different random errors E. 

Note that the simulations of simple and improved decoding differ in how the third 
step is executed. In particular, in the third step, simple decoding guesses an X error in 
the north-west corner of each square with non-trivial syndrome D, whereas improved 
decoding uses also syndrome B to guess the specific location of the X error within 
the square with non-trivial D syndrome, as described in Section 14.31 Note that both 
decoders correct Z errors at location 4 in the unit cell (see Figure 10(b) ) in the fourth 
step if the corresponding syndrome B is non-trivial. 

The minimum weight matching in the fifth step of the algorithm is used as follows. 
Since the fifth step is repeated twice, once using syndromes depicted in Figure [TTl and 
then analogously using syndromes depicted in Figure [I2l we only describe the former. 
We identify the syndromes A and C in Figure [TT] that are non-trivial, and represent each 
such syndrome as a vertex. Then we calculate the pairwise distance of all such vertices 
and use the minimum weight matching algorithm to pair up vertices so that the sum of 
the distances between the vertex pairs is minimized. The distance of syndrome u and 
V in the graph represents the minimal number of Z errors that have to be introduced 
to the original lattice in Figure [TT] to change syndromes u and v while keeping all other 
syndromes intact. The Z errors that need to be introduced to change syndromes u and 
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Figure 15. The threshold is around 1.5%. 
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Figure 16. The threshold is around 2%. 

V lie on the shortest hne connecting the two syndromes. These Z errors are added to 
Eguessed the fifth step of the simulation. 

We executed the experiment with the depolarizing error model with error 
probability ranging from 0% to 5% with 0.1% increments. The experiment was repeated 
for square lattices consisting of a total of 640, 2, 560, 10, 240, and 40, 960 qubits. Each 
of these square lattices consisted of n by 2n unit cells where each unit cell has 20 qubits. 
For example, the lattice in Figure [TT] has 2 by 4 unit cells for a total of 160 qubits. For 
each considered depolarizing error rate and lattice size, we repeated the error correction 
1, 000 times and recorded the number of experiments that resulted in a failed decoding. 

The results are depicted in Figure [15] for the simple decoding algorithm, and 
Figure [16] for the improved decoding algorithm. We observe that in both cases the 
code exhibits a threshold. In the former case the threshold is around 1.5%, and in the 
later case around 2%. 

As the error rate increases, the failure probability in Figure [I5] and [16] approaches 
15/16, which is as expected since the success probability for picking the correct coset at 
random is 1/16. 

We did not perform experiments with the square-octagon code as this code and 
its noise threshold with noiseless error correction are studied in [31j. Surprisingly the 
authors find a similar depolarizing noise threshold of 2%. 
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6. Conclusion 

We studied tlie error correction properties of topological subsystem codes whose gauge 
group has local two-qubit generators. We provided an example of such code on a 
two-dimensional grid of qubits and we described a candidate decoding algorithm. Our 
Monte Carlo simulations indicated that our code exhibits a threshold for a depolarizing 
error rate of around 2%. This is substantially lower than the threshold for the surface 
code with noise-free error correction. We have tried to assess whether the threshold 
achieved by our improved decoder is close to the theoretical, achievable, threshold in 
Eq. ( 1T6|) . For this we have numerically executed a simulated annealing algorithm that 
performs minimum-weight decoding. The algorithm solves the problem of finding the 
minimum energy of the 2D local Hamiltonian with quenched disorder determined by 
PE by gradually lowering the temperature. It then compares the found minima for 
various P and picks the one with the lowest minimum. Our numerics indicate that 
this decoder does not improve on the decoders described in the previous section. It 
is an important open question whether there are other topological subsystem codes 
which have a threshold which is closer to that of the surface code or whether the lower 
threshold is an inevitable price that we pay by enhancing the locality, 2-body versus 
4-body, of the code. 
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Appendix A 

In this section we prove that an eigenvalue of a stabilizer S can be determined 
using eigenvalue measurements of the gauge group generators Kj only if S admits a 
decomposition Eq. ([5]) satisfying Eq. ([6]). 

Indeed, the most general measurement protocol involves adaptive eigenvalue 
measurements: suppose we have already measured eigenvalues of some generators 
Ki, . . . ,Kj obtaining eigenvalues 71, . . . , 7^- G {+1,-1}. The next generator to be 
measured, Kj^i, may be chosen as an arbitrary function of the outcomes 71,..., 7^ 
and, may be, some extra random bits. Let 
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be the projector describing the eigenvalue measurement of Kj. Let 
Rj = Uj---Ui 

be the operator describing the entire sequence of measurements up to the step j for 
some fixed configuration of outcomes 7 = (71, . . . ,7m)- We shall say that the above 
adaptive measurement protocol simulates eigenvalue measurement of S if there exist a 
function a : {+1, -1}™ {+1, -1} such that 

SR^ = R^S = a{-i)R^ (17) 

for all 7 G {+1,-1}™. Equivalently, the reduced state obtained after the last 
measurement must be an eigenvector of S with an eigenvalue (7(7) . Below we prove 
that this is possible only if S admits a decomposition Eq. ([5]) satisfying Eq. (|6]), may 
be, with a different m. 

Indeed, for any fixed outcomes 7 define a family of stabilizer codes (Abelian groups 
of Pauli operators) Sq, Si, . . . , Sm C Vn such that Sq = (/), Si = {■yiKi), and 

= (7,+iir,+i, 5;), (18) 

where 

S'^ = {Pe S, : PK,+i = K,+iP} (19) 

for j = 1, . . . , m — 1 and = Sm- Using the standard stabilizer formalism [33] one 
can easily check that Rj = Iljf/j, where Ilj is the projector onto the codespace of 
the stabilizer code Sj and Uj is some Clifford group unitary operator that belongs to 
the algebra generated by Ki, . . . , Kj. In particular, both 11^ and Uj commute with 
S. Obviously, condition Eq. ( ITTl) holds only if a{'y)S G Sm- Using the definitions 
Eqs. ffTMOD we conclude that 

cT{^)S={^mKmr^^^Sm-i for somc Sm-ieS'^,, (20) 

and some a(m) G {0, 1}. Note that Km commutes with Sm-i by definition of the group 
S'm-i- Using the inclusion Sj C Sj inductively we arrive at the decomposition 

a(7)5 = i^mKmr^""^ ■ ■ ■ iliKiT^'^ (21) 

for some a(l), . . . , a{m) G {0, 1}. This decomposition obeys {'jjKj)"'^^'' ■ ■ ■ {'jiKi)°'^^^ G 
Sj for all j = 1, . . . , m. Hence K°'^''^ commutes with i^'"!^"^^ ■ ■ ■ Ki^^'' for all j. Omitting 
all the factors with a(j) = we arrive at the decomposition Eqs. ( !5]l6l) . 
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